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Abstract 

An experimental study was conducted on two small rockets (1 ION thrust class) to directly compare a standard conical nozzle with 
a bell nozzle optimized for maximum thrust using the Rao method. In large rockets, with throat Reynolds numbers of greater than 
lxlO 5 , bell nozzles outperform conical nozzles. In rockets with throat Reynolds numbers below lxl 0 5 , however, test results have 
been ambiguous. An experimental program was conducted to test two small nozzles at two different fuel film cooling percentages 
and three different chamber pressures. Test results showed that for the throat Reynolds number range from 2x1 0 4 to 4x1 0 4 , the bell 
nozzle outperformed the conical nozzle. Thrust coefficients for the bell nozzle were approximately 4 to 1 2 percent higher than those 
obtained with the conical nozzle. As expected, testing showed that lowering the fuel film cooling increased performance for both 
nozzle types. 


Introduction 

Future space exploration will require increasing payload 
fraction and decreasing system masses. Therefore, maximizing 
rocket engine performance is an important design goal. Maxi- 
mizing the thrust for a rocket engine can be achieved by 
optimizing the nozzle contour. Optimizing the performance of 
finite length nozzles is typically accomplished using an invis- 
cid core flow and a boundary layer displacement. G.V.R. Rao 
developed a method which optimizes a rocket nozzle contour 
for a given length or expansion ratio such that maximum thrust 
is achieved. 1 Rao’s method was based on the assumption of 
inviscid isentropic flow. 

This method has been used in many studies of different classes 
of rocket engines with a variety of results. Rocket engines are 
generally classified by their throat Reynolds number based on 
uniform flow and properties at the throat. For engines with 
Reynolds numbers larger than lxlO 5 , the Rao optimization is 
quite effective, because the ratio of boundary layer flow to the 
total flow is small. In a study by Farley, 2 three large Rao 
optimized bell nozzles were compared to a 15 degree conical 
nozzle. The thrust produced by the optimized nozzles was greater 
than that obtained with the conical nozzle. In fact, bell contour 


nozzles have been used routinely for many years in large liquid 
rocket engines. 3 ’ 4 ’ 5 Conical nozzles are typically used only when 
fabrication and design costs outweighed performance. 

In very small thrusters with Reynolds numbers less than 
lxlO 4 , viscous boundary layer effects are large in comparison 
to the total flow and a bell contour is not nearly as effective, as 
discussed in a pair of experimental studies. 6 ’ 7 In these studies, 
a variety of small nozzles were tested to determine the effect of 
contour, propellant, expansion ratio, and Reynolds number on 
performance. Normally, factors in performance loss are diver- 
gence, boundary layer effects, and heat transfer effects to the 
walls. By using ambient temperature gas, the heat transfer 
effect was eliminated in both studies. The nozzles were conical, 
trumpet, and Rao optimized bell contours. One paper con- 
cluded that the difference in contours had no effect on either 
divergence or viscous losses. 6 The other concluded the bell 
contour had the lowest performance of the nozzles tested. 7 The 
reasons given in both were that the viscous effects of boundary 
layer were so large that an inviscid isentropic assumptions were 
not valid. 

For rockets that are sized between the very small and large 
ones the usefullnes of the Rao optimization code is less clear. 
Investigations on a gaseous hydrogen/oxygen engine were 
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performed by M. A. Appel with different nozzle contours in the 
134 to 445 N thrust class. 8 A rhenium 80 percent bell nozzle, a 
regeneratively cooled nozzle, and a stainless steel 1 8 degree half 
angle conical nozzle were each tested using the same injector in 
this effort. The stainless steel nozzle was built to verify injector 
and thrust stand operation prior to testing the other two nozzles. 
The regeneratively cooled nozzle had different dimensions for a 
different aspect of the study, but the conical and bell nozzles had 
the same length, exit area and expansion ratio. Unexpectedly, the 
test results showed that the conical nozzle provided a higher 
specific impulse than the bell nozzle. The author postulated that 
the conical nozzle operated with a higher kinetic efficiency and 
suggested this as an area for further investigation. Kushida 9 
examined a 22 N monopropellant hydrazine thruster with a 1 00: 1 
area ratio conical nozzle, a 300: 1 Rao optimized bell, and the same 
bell nozzle cut off at expansion ratios of 200: 1 and 1 00: 1 . Results 
showed that the 300:1 bell nozzle outperformed the conical 
chamber but the conical was the better performer than the lower 
area ratio bell nozzles. A direct comparison was difficult, how- 
ever, because several dimensions were changed in addition to the 
nozzle contours. Both the length and expansion ratio were greater 
for the optimized nozzle, which generally improves performance 
in a rocket engine. Based on the ambiguous results reported to 
date, the validity of the Rao nozzle optimization method is 
uncertain for thrusters with Reynolds numbers between lxlO 4 
and lxlO 5 . 

In the study reported herein, testing was performed to di- 
rectly evaluate the applicability of Rao optimization methods to 
rockets with Reynolds numbers below lxlO 5 . The performance 
of two 1 10 N chamber/nozzles with throat Reynolds numbers 
between 2xl0 4 and 4xl0 4 was measured to examine the effect 
of nozzle contour on performance in this intermediate Reynolds 
number range. Both nozzles were designed with the same 
chamber profile, throat diameter, and expansion ratio. Because 
the chamber and nozzle are one continuous structure, each 
piece of hardware is simply identified by the nozzle shape, ie. 
bell or conical, from here on. The conical nozzle had a 19 degree 
expansion half angle. This rocket was originally tested in 
1993. 10 The bell nozzle was designed via the Rao code. As part 
of the effort, a one dimensional equilibium code, ODE, was 
used to generate theoretical performance predictions to calcu- 
late the characteristic velocity efficiencies presented later in 
this study. 11 

Test Hardware 

Both nozzles tested were fabricated via the same techniques 
using oxygen free high conductivity (OFHC) copper with water 
cooling passages in the walls. The inner housings were milled 
first to cut water cooling passages into the outer surface. The 
chamber profiles were identical from the injector face to the 
throat. The contours of the nozzles are shown in figures 1 (a) and 
(b), and further details of the bell nozzle contour are given in 


table I. The divergence angle of the conical nozzle was about 
19 degrees and the divergence angle of the bell nozzle was 
about 9 degrees at the exit plane. The design specified throat 
area was 1.267 cm 2 . Post fabrication measurements of the 
chamber, throat, and exit gave dimensions that were very 
similar as shown in table II. The resulting dimensional differ- 
ences based on area were 0. 10 percent at the chamber, 2.2 per- 
cent at the throat, and 0.14 percent at the nozzle exit. 

The outer housing of each nozzle was machined in two pieces 
separated along the thruster axis. The sections fit around the 
inner housing as a shell to insure the interior contour was one 
continuous profile from chamber inlet to nozzle exit. The outer 
housing halves were attached to each other with bolts around 
the inner housing. All seams were filled with high temperature 
solder. Both chambers were instrumented with thermocouples 
and pressure taps in several quadrants along the contour. 

The same injector, 12 shown schematically in figure 2, was 
used in all testing. The injector was designed to produce a core 
flow of combustion products with an annulus of fuel film 
cooling (FFC). In the core, oxygen was injected through a 
platelet stack into the combustion chamber just upstream of a 
spark plug used for ignition. The hydrogen flow was divided by 
means of a fuel splitting washer, with part injected downstream 
of the spark plug into the core and the rest used to fuel film cool 
the chamber wall. The washer was specifically designed to give 
a predetermined fuel film cooling split from the total fuel flow. 
The film cooling was directed along the chamber walls by an 
annular sleeve. By changing the washer, the amount of film 
cooling could be varied between 55 to 85 percent. A 61 percent 
FFC percentage was selected as the baseline for this study since 
that was nominal design operating point for this injector. Tests 
at 75 percent FFC were also run to investigate the sensitivity to 
FFC changes. 

Test Facility 

All tests were conducted in a small rocket test facility, 13 in 
which an equivalent vacuum of 36.6 km altitude is achieved 
using air ejectors. Gaseous hydrogen and gaseous oxygen were 
used for propellants. A schematic is shown in figure 3. The tank 
is 1 .8 m long between the flanges and 0.9 m in diameter. Inside 
the tank, the rockets are oriented horizontally to fire through a 
diffuser into an exhaust quenching chamber. The rocket 
exhaust is pulled through the pair of ejectors and vented to 
atmosphere by mufflers. The ejectors are supplied with 
900 kPa air at 4.5 kg/sec to provide the pumping for the altitude 
tank. 

Thrust is measured by mounting the rocket on a thrust stand. 
The stand and rocket are supported with flexure plates. Both 
float free on the flexures and apply a load to a strain gage load 
cell through a ball joint. To ensure accurate thrust calibration, 
the ball joint must make contact with the load cell at all times. 
The thruster is positioned to apply a small load on the load cell 
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even when the thruster is not firing, to ensure contact with the 
load cell is maintained. The small load is accounted for in thrust 
calibration and is subtracted from the rocket thrust. The propel- 
lant and pressure transducers lines are mounted perpendic- 
ular to the rocket axis with stainless steel tubing to provide 
repeatable thrust measurements. The thrust stand is calibrated 
by applying loads through a piston assembly to compare the 
thrust measurement load cell against another previously cali- 
brated reference load cell. Pretest, posttest, and applied load 
tests are taken at altitude conditions with pressurized propellant 
lines. All thrust measurements are corrected for background 
pressure to determine vacuum thrust. 

Rao Optimization Code 

Rao discussed the need for nozzle contour optimization. 1 4 In 
order to maximize thrust, a thruster is designed to produce 
uniform parallel exit flow, with the exit plane pressure equal to 
the ambient pressure. In a conical nozzle, flow divergence 
losses are significant, particularly as cone angles increase. Rao 
suggested that contouring the nozzle wall would turn the flow 
closer to the thruster axis and thus reduce flow divergence. He 
developed an optimization method for designing nozzles which 
assumes inviscid isentropic flow. 15 The method uses the calcu- 
lus of variations for an ideal gas with constant gamma expan- 
sion to calculate the optimum nozzle contour. Supersonic flow 
properties across the nozzle throat are input to start the solution. 
Transonic flow in the throat region is calculated as a function 
of the upstream radius of curvature at the throat and the specific 
heat ratio for axisymmetric flow. The method iterates to an 
optimum nozzle contour defined to meet for one of three initial 
conditions: a given length, a given expansion ratio, or a fixed 
envelope (exit radius and length). The Rao code was used in this 
study to obtained an optimized bell nozzle contour with the 
dimensions shown in table I using the fixed envelope condition. 

Test Procedures 

The same tests were conducted on each nozzle using gaseous 
hydrogen and oxygen propellants. Tests were conducted at 
three different chamber pressures, 500, 370, 255 kPa, and two 
FFC percentages, 61 and 75 percent, over a mixture ratio range 
of 4 to 8. The design condition for the rocket hardware was 
500 kPa with 61 percent FFC. The higher FFC percentage and 
the lower chamber pressures were chosen to allow testing 
through a throat Reynolds number range of 2x 1 0 4 to 4x 1 0 4 . The 
throat Reynolds number (Re) range was calculated using the 
equation 


where m is the mass flow rate of the propellants, p is the 
absolute viscosity of the combustion products, and D t is the 
diameter of the nozzle throat. The viscosity is determined from 
theoretical predictions generated by the Gordon-McBride 
program. 11 The two percent difference in throat areas some- 
times resulted in slightly lower chamber pressure in the conical 
nozzle. A typical test began by pumping the test tank to altitude 
and pressurizing the propellant lines with fuel. A thrust stand 
calibration was then performed. After the calibration, a con- 
stant load was applied to the thrust stand, which provided data 
for determining the uncertainties. Once the applied load was 
removed and inlet pressure conditions at the injector were 
adjusted, the rocket engine was fired. By modifying the propel- 
lant inlet pressures, mixture ratios from 4 to 8 were tested. Once 
tests were completed for the day, another constant load test and 
another load cell calibration were conducted. 

Measurement Uncertainties 

Measurement uncertainties were calculated using the 
JANNAF recommended procedure. 16 The uncertainties were 
calculated from the precision (random) and bias (constant) 
errors associated with the measurement instruments and data 
acquisition system. Precision errors propagated through the 
performance parameters of specific impulse, characteristic 
velocity, and thrust coefficient are shown in table III. Initial 
testing with the bell nozzle for the 500 and 370 kPa chamber 
pressure cases at 61 percent FFC had a very small pretest load 
(less than 0.5 N) on the thrust stand prior to ignition. Since the 
thrust stand could float free of the load cell at this load, a larger 
load was applied prior to ignition to ensure pretest contact 
between the thruster and the load cell. The low pretest load 
affected the repeatability of the thrust stand zero. In an effort to 
ensure that the thrust measurement was not affected by the 
small pretest load, the pretest load was increased and a few 
selected points in that test series were repeated. A field check 
at that time showed performance results consistent with the 
original tests at the smaller pretest load. However, in subse- 
quent data reductions, the tests with the small pretest load had 
significantly larger uncertainties. The range of uncertainties for 
both cases are shown in table HI, the cases with the smaller 
pretest load being indicated by parenthesis. The remainder of 
the test data uncertainties were also presented in the same table. 
In general, the greatest uncertainties were at the lower chamber 
pressures. 

Results and Discussion 

Injector/Chamber Performance 

To show that the injector and chamber behaved identically 
for the two nozzles, the characteristic velocity efficiencies are 
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plotted in figures 4 and 5 as a function of mixture ratio for tests 
run at 61 and 75 percent FFC, respectively. The efficiencies 
were determined by dividing the experimental characteristic 
velocity by the theoretical equilibrium characteristic velocity 
predictions. 11 As expected, the plots show no measurable 
differences in this parameter between the bell and conical 
nozzles as the injector and chamber profiles are the same. The 
slight difference in throat area resulted in slightly higher 
chamber pressures (approximately 1 percent), in the bell nozzle 
for the same inlet conditions. Clearly the difference did not 
have a measurable impact on the performance of the injector 
and chamber, since the characteristic velocities measured with 
both nozzles fell within the same uncertainty band. The overall 
efficiency is slightly higher for the 61 percent FFC cases. This 
was an expected result because more hydrogen is reacted in the 
core combustion flow at this condition. Similiar trends are 
shown in the plots of characteristic velocity versus mixture 
ratio (figs. 6 and 7). 

Nozzle Performance 

Vacuum specific impulse versus mixture ratio and thrust 
coefficient versus mixture ratio are typically used to illustrate 
nozzle performance. The vacuum specific impulse versus mix- 
ture ratio for the 61 and 75 percent FFC cases are plotted in 
figures 8 and 9, respectively, for each of the three chamber 
pressures tested. The data were fit to a second order polynomial 
curve and the measurement uncertainties added. At 6 1 percent 
FFC (fig. 8), the plots of nozzle data show that the bell nozzle 
specific impulse are at least 5 percent higher at 500 kPa, 
6 percent higher at 370 kPa, and 12 percent higher at 255 kPa 
than that of the conical for the 61 percent FFC cases. In 
75 percent FFC cases (fig. 9) the bell nozzle specific impulse is 
at least 5 percent higher at 500 kPa, 10 percent higher at 
370 kPa, and 14 percent higher at 255 kPa. These results are 
statistically significant and show clearly that the bell nozzle out 
performed the conical nozzle. The higher performance obtained 
with lower FFC percentage was expected as discussed above. 

Thrust coefficients are plotted for the 6 1 and 75 percent FFC 
cases in figures 1 0 and 1 1 , respectively, again for each chamber 
pressure case. As seen in the specific impulse plots, the bell 
nozzle had the higher performance for both FFC percentages. 
The thrust coefficients at 61 percent FFC are generally 6 per- 
cent higher at 500 kPa, 7 percent higher at 370 kPa, and 
1 1 percent higher at 255 kPa for the bell nozzle than for the 
conical. The thrust coefficients at 75 percent FFC are generally 
4 percent higher at 500 kPa, 8 percent higher at 370 kPa, and 
1 1 percent higher at 255 kPa for the bell nozzle than for the 
conical. Again clearly, the bell nozzle outperforms the conical 
nozzle for this Reynolds number range of between 20,000 to 
40,000. 


Conclusions 

Testing was conducted to directly compare the performance 
of a bell nozzle designed by the Rao optimization method and 
a conical nozzle of the same size in the 2X10 4 to 4X10 4 throat 
Reynolds number range. Previous studies had indicated that an 
optimized bell nozzle may not outperform a conical nozzle at 
these Reynolds numbers. To eliminate the effects of injector 
performance, the same injector and chamber profile were used. 
Tests were conducted over a mixture ratio range of 4 to 8 for 
three different chamber pressures and two fuel film cooling 
levels. Characteristic velocity measurements indicated that the 
injector and chamber gave the same performance with both 
nozzles. Plots of vacuum specific impulse and thrust coeffi- 
cients indicated a higher performance by the bell nozzle. Also, 
higher performance at the lower FFC percentages for both 
nozzles. For the throat Reynolds number range of this design, 
the test results confirm that the Rao optimized bell nozzle yields 
higher performance than a simple conical nozzle design. 
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TABLE 1— WALL CONTOUR POINTS FOR 
RAO NOZZLE BEGINNING AT THROAT 
TO NOZZLE EXIT 


Axial 

distance, 

cm 

Radial 

distance, 

cm 

Axial 

distance, 

cm 

Radial 

distance, 

cm 

0.000 

0.638 

1.267 

1.349 

0.157 

0.681 

1.359 

1.400 

0.175 

0.691 

1.450 

1.450 

0.203 

0.709 

1.544 

1.501 

0.231 

0.726 

1.623 

1.542 

0.262 

0.744 

1.727 

1.595 

0.295 

0.765 

1.910 

1.687 


0.785 

2.217 

1.836 

0.366 

0.808 

2.621 

2.017 

0.404 

0.833 

2.891 

2.129 

0.444 

0.856 

3.152 

2.233 

0.488 

0.884 

3.495 

2.362 

0.533 

0.912 

3.843 

2.487 

0.579 

0.942 

4.280 

21631 

0.632 

0.973 

4.740 

2.774 

0.686 

1.006 

5.314 

21936 

0.744 

1.041 

5.928 

3.096 

0.805 

1.080 

6.302 

3.185 

0.871 

1.120 

6.891 

3.315 

0.940 

1.161 

7.534 

3.444 

1.016 

1.204 

8.037 

3.538 

1.095 

1.252 

8.434 

3.607 

1.179 

1.298 

8.989 

3.693 


TABLE II — CROSS SECTIONAL AREAS OF THE 
BELL AND CONICAL NOZZLES 


Location 

Cross sectional area, 
cm 2 , cm 

Bell nozzle, 
radius 

Conical nozzle, 
radius 

Chamber 

5.045 (1.268) 

5.050 (1.267) 

Throat 

1.248 (0.637) 

1.275 (0.630) 

Nozzle exit 

42.65 (3.687) 

42.70 (3.685) 


TABLE III.— PERFORMANCE PARAMETER UNCERTAINTY PERCENTAGES 
FOR THE BELL AND CONICAL NOZZLES AS A FUNCTION 

OF CHAMBER PRESSURE 


Uncertainty, ± percent 

Bell nozzle 

Conical nozzle | 



61% FFC 

75% FFC 

61% FFC 

75% FFC 

Specific 

500 kPa 

16 

16 

1.1 

1.1 

impulse 

370 kPa 

(4.0) 

3.3 

19 

1.4 

1.5 


255 kPa 

(6.8) 

5.7 

4.6 

10 

1.9 

Characteristic 

500 kPa 

1 2 

1.1 

1.2 

1.1 

velocity 

370 kPa 

1.6 

1.4 

1.5 

1.5 


255 kPa 

2.4 

2.1 

10 

10 

Thrust 

500 kPa 

2.6 

2.6 

0.4 

0.5 

coefficient 

370 kPa 

(4.1) 

3.3 

2.9 

0.7 

0.7 


255 kPa 


4.5 

1.1 

0.9 
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Figure 4. — Characteristic velocity efficiency versus mixture ratio for 
61 percent FFC for both bell and conical nozzles. 
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Figure 5. — Characteristic velocity efficiency versus mixture ratio for 
75 percent FFC for both bell and conical nozzles. 
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Figure 6. — Characteristic velocity versus mixture ratio for 61 percent 
FFC for both bell and conical nozzles. 
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Figure 7. — Characteristic velocity versus mixture ratio for 75 percent 
FFC for both bell and conical nozzles. 
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Figure 8. — Specific impulse versus mixture ratio for 61 percent FFC for 
both bell and conical nozzles, with uncertainty bars. 
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Figure 9. — Specific impulse versus mixture ratio for 75 percent FFC for 
both bell and conical nozzles, with uncertainty bars. 
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Figure 1 0. — Thrust coefficient versus mixture ratio for 61 percent FFC for 
both bell and conical nozzles, with uncertainty bars. 
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Figure 1 1 .—Thrust coefficient versus mixture ratio for 75 percent FFC for 
both bell and conical nozzles, with uncertainty bars. 
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